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ABSTRACT 

Particles inside a planetary ring are subject to forcing due to the central 
planet, moons in inclined orbits, self-gravity of the ring and other forces due 
to radiation drag, collisional effects and Lorentz force due to magnetic field 
of the planet. We write down the equations of motion of a single particle 
inside the ring and solve them analytically. We find that the importance of 
the shear caused by variation of the radial velocity component with local 
vertical direction cannot be ignored and it may be responsible for damping 
of the bending waves in planetary rings as observed by the Voyager data. We 
present the wave profile resulting from the dissipation. We estimate that the 
surface mass density of the C ring to be of the order of a ~ 1.2 — 1.6gm cm -2 , 
and the height h ~ 2.2 — 2.4m. These theoretical results are in agreement with 
observations. 
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1 INTRODUCTION 

Voyager photo-polarimeter (PPS) experiment indicated that the thickness of the Saturn's 
ring is probably smaller than 200m at the outer edge of the A -ring (Lane et al, 1982). 



* E-mail: b_abhijit2k@yahoo.co.in 
f E-mail: chakraba@bose.res.in 



2 Abhijit Bhattacharyya and Sandip K Chakrabarti 

When the oscillatory features in bending waves are identified with possible resonance events 
(Shu et al, 1983), after a fit of the ring profile, the effective thickness is found. The effective 
thickness of the ri ng is found to be only a few tens of meters (Rosen et al. 1988, 1991), 
Chakrabarti (1989) and Chakrabarti and Bhattacharyya ((20011 1 . One of the important new 
ingredients that has gone into estimation of the thickness is a new source of shear which 
was found to be present during detailed numerical study of dynamics of particles inside the 
Saturn's ring. This extra source of shear is found to effectively damp out the 5:3 bending 
wave of Mimas within a few tens of kilometers and since then it is considered to be important 
(Nicholson et al 1991, Rosen et al, 1991 and Brophy et al, 1992) to describe correctly the 
particle dynamics inside a ring. 

We now present briefly so me of the work done on vertical structure and wave profile prior 
to ours. Simon and Jenkins (1994) performed their analysis assuming the system having a 
dilute flow maintained by identical smooth spherical frost-free particles. The balance laws 
were obtained by taking moments of the Boltzmann equation. While taking care of collisions 
between particles, they assumed the normal velocity components, before and after the col- 
lision, to be related through the coefficient of restitution, while the tangential components 
were unaffected. They also assumed a symmetry about the equatorial plane (z = 0) and that 
there was no mean motion normal to the mid-plane. They concentrated on the outer edges 
of the Saturn's A -ring and used parameters near #:7Lindblad resonance to obtain relations 
between granular parameters and the co-efficient of restitution and optical depth. They find 
the coefficient of restitution to be monotonically increasing with the optical depth. No effort 
was mad e to g ive the profile of the bending wave. 

Salo ()1991l ) studied viscous stability properties of dense planetary rings with numeri- 
cal simulations. He confirmed earlier results of Wisdom and Tremaine (jl988T ) that for the 
standard elastic model of icy particles, the viscous instability is not expected for identical 
meter-sized particles, while for denser systems with centimeter sized particles the possibility 
persists. Salo also did simulations including self-gravity resulting in the variation of particle 
number density with vertical height. In this simulation, Salo considered 100-200 impacts 
per particle. He assumed that the particle density is the highest on the mid-plane of the 
ring. This simulati on ma y be applicable for Saturn's E, F and G ring, Uranian and Jovian 
rings. Later, Salo |l992l ) did numerical simulations of collisional systems with power law 
distribution of particle sizes and found that equilibrium geometric thickness of the Saturn's 
ring to be around 25m for layers of cm-sized particles, and around 10m for larger particles. 
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Schmidt et al. ()1999l ) investigated viscous oscillatory over- stability of an unperturbed 
dense planetary ring which may play a role in the formation of ra dial structure in Saturn's 
B ring. They extended the modeling of Simon and J enkins (|l9941 ) and compared Ga ussian 
ansatzs which were used in Goldreich & Tremaine (|l978l ) and Araki & Tremaine ((1986) 
with hydrodynamical ansatzs. They considered enhanced vertical frequency which appears 
to flatten the disk and enhances collision frequency. The simulation was carried out with 
identical sized particles. This simulation shows initial onset of over-stability in the ring and 
leads to the development of radial structure although no explanation was presented on the 
damping of the bending wave. 

In our paper, we are considering approximate analytical equations obtained from the 
kinematic formulation for Titan -1:0 vertical resonance to estimate a relation between surface 
mass density, ring height and the damping length of the bending wave in the C -ring. Our 
approach emphasizes the effect of vertical motion of particles over and above the epicyclic 
one. While doing the vertical and horizontal excursion, the particles collide and transport 
vertical component of momentum, thereby damping the vertical profile of the wave. In the 
present paper, we write down the equations governing the particle dynamics including this 
form of shear. We compute analytically the shear developed for Titan -1:0 bending resonance 
and then study the damping length as a function of surface mass density and height as free 
parameters. We restrict our attention to that region of the parameter space where the 
properties of the wave profile is comparable to those observed. There are degeneracies in the 
parameter space in the sense that the same damping length could be achieved for high and 
low surface mass densities (a). Upon scrutiny it was observed that one of the solutions (that 
corresponding to lower a) created too many wave profiles than what is observed and thus 
it was rejected. The present analysis provides the wave profile and therefore a comparison 
with the observation is easier to make. Our earlier (Chakrabarti and Bhattacharyya, 2001) 
procedure could give population distribution throughout the ring at different phases which is 
not possible in the present method. A realistic three dimensional simulation would probably 
be essential to understand the correct description of the problem. 

It is to be noted that we have included only the effect of energy loss due to infrequent 
collision among shearing layers and no attempt has been made to include cooling effects 
because of poorly understood physics. The collisional damping reduces kinetic energy and 
the degree of vertical excursion. This is used at each radius to obtain the wave profile. We 
computed the coefficient of restitution in our ring and found this to be close to unity, in near 
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Figure 1. Choice of the co-ordinates in our calculations. The broad view with the planet and the ring is shown in (a). In 
(b) we show details of the co-ordinates used in the small box of (a). Here, the mid-plane of the bending wave is shown to be 
flanked by two artificial 'boundaries' located at ±h distance from the mid-plane. The particle P is located at (x',y',z') in the 
co-ordinate system with origin at (x*, y*, z*). See text for details. 



agreement with the results of Simon and Jenkins (|1994l ) where the coefficient is found to be 
0.9 for the same disk parameter. In other words, we assumed nearly-elastic collision among 
the particles. In future, we shall model the physics of breakup and merger of 'particles' in 
the ring and other associate effects. 

In the next section, we present the basic equations of the motion of particle and explain 
the genesis of each of the terms. In §3, we present the solutions of the equations and compute 
the average shear analytically. In §4, we apply our results in various ring conditions. Finally, 
in 55, we draw our conclusions. 



2 MODEL EQUATIONS 

We choose the usual right handed Cartesian coordinate system (X, Y, Z) at a radial distance 
r from the center of the planet with the origin located on the equatorial plane of the planet 
with X — Y plane coinciding with this plane (See, Fig. Q). 

The X -axis points radially outward and Y -axis points toward the azimuthal direction. 
The frame is rotating around the planet with the local Keplerian frequency fl(r). Let the 
amplitude of the bending wave be e. The mid-plane itself oscillates up and down with this 
amplitude. Matter oscillates up and down around the instantaneous mid-plane with an 
amplitude h! — h/2 - half of the thickness h of the ring. Let the coordinate of the origin O' 
of a Cartesian frame (X', Y', Z') which is oscillating with the mid-plane of the disk be (x*, 
y*, z*). In the absence of oscillations of the mid-plane, these two coordinate systems merge. 
A particle moving within the ring having coordinate (x, y, z) has a coordinate of x\=x-x*, 
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yi=y-y* and zi=z-z* with respect to O' . If uj be the disturbance frequency due to moon 
forcing, i.e., the angular frequency of the propagating wave, then the phase of the mid-plane 
is 4>*=k x x*+k y y*-ut and that of the particle located at a point P (x, y, z) is 4>=k x x + k y y - 
ujt (Here k x and k y are the x— and y— components of the wave vector k.) so that the phase 
difference between the particle and the midplane is: 

(j) - (f>* = k x (x - x*) + k y (y- y*) = k x x x + k y y 1 . (1) 

The frequency u has a linear relation with orbital frequency of moon around the planet 
(Qm), vertical (hm) and epicyclic (%) frequencies of moon given by: 

uj = mVL M ± n^ M ± P^m, ■ (2) 

where, m, n and p are non-negative integers. For small perturbations from a satellite, the test 
particle responds as a multi-dimensional, forced, linear-harmonic oscillator. In particular, it 
may suffer resonances if the relative frequency at which it experiences the disturbance of 
the satellite referred to its local rotation rate, is equal to a ny of the natural frequencies 
of its free oscillations as discussed by Franklin and Colombo (1970). The test particle will 
suffer vertical resonance due to inclination of orbit of satellite if its radial distance of vertical 
resonance (ry) from the center of the planet satisfies 

uo - mfi(ry) = ±/i(r y ). (3) 

We can always consider Qm < because the large satellites orbit outside the main ring 
system. Here, m = Q(r)/Q M designates a particular resonance (See, e.g., Shu, 1984 and 
Murray and Dermott, 1999 for details.). Satellites will impart forcing of different strength 
at different orbits and this is taken care of by putting different values of m, n and p in 
Eq. El We considered Titan -1:0 resonance for which m — 1 in used in Eq. El This mode 
is the strongest resonance imparted by Titan on Saturn's ring. By definition, uj = (flp-fl), 
where, Qp is the pattern frequency. In the present resonance, this is the angular frequency 
of the perturbing moon with an opposite sign of orbital frequency of moon. It is to be kept 
in mind that Titan has a semi-major axis of about 1,221,850 km to revolve round Saturn 
along an orbit having an orbital inclination of 0.33 degree with Saturnian equatorial plane. 
Due to this inclination of orbit, spiral bending wave is launched in the ring at the resonance 
location ry = 77514.8 km which is inside Saturn's C ring. Hence the properties of the C 
ring could be obtained by understanding the behaviour of the bending wave due to Titan 
-1:0 resonance. 
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Let (x', y', z') be the coordinate of the particle in the (X', Y', Z') frame. Thus, 

z' = h + ecos0*, (4) 

where, h is the half-thickness of the ring as stated earlier. It is easy to show that the local 
normal at (x 1 , y', z') on the mid-plane is given by, 

(5) 



n = 



; ek x sin <p* 



J 



ek y sin (j>* 



V 



w J w 

where, W^=(l + e 2 k 2 sin 0*) 1//2 . From this, one derives the following important relations, 

xi = z\ek x sin0*, 
and y 1 = Ziek y sm(p*. 



(6) 
(7) 



It is conventional to write the potential (f>p due to the planet in the equatorial plane as 



(r,0) = 



GMf 



(8) 



l-J2J2n(R P /ry n P 2n (0) , 

n=l 

where, M P is mass of the planet, R P is its radius, J 2n is its 2n th multipole moment, P2n(0) 
is the Legendre polynomial of order 2n at the origin. If one ignores the oblateness of the 
planet, then the multiple moment terms are absent and the vertical frequency fi of the 
particle defined by, 

2 - 8 ^ P (9) 



// 



dz 2 



becomes identical to the Keplerian frequency. The epicyclic frequency k of the particle 
defined by, 



K 



1 d 

r 3 dr 



(r 2 n) [ 



(10) 



This is also identical to the local Keplerian frequency Q. 
In general, there will be three major acceleration terms, 

9 = 9p + 9d + 9m, (11) 

where, subscripts P, D and M denote the acceleration due to the planet, the self-gravity 
of the ring and the moon respectively. In the first approximation, one can assume that the 
vertical motion is due to the planet and the moon only, so that the vertical component of 
the equation of the motion is given by, 

z = —Q 2 z + a cos (k x x + k y y — out) . (12) 
The solution of this equation is, 

n 

cos (k x x + k y y — cot) . (13) 



z = 



tt 2 -tu 2 
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The factor in front of the cosine term can be identified with the amplitude e of the vertical 
movement of the mid-plane due to the satellite forcing provided e (Q 2 — uo 2 ) = a. Thus the 
forcing term due to the satellite is, 

(14) 



Qm = e (p — 00 j cos (k x x + k y y — cut) . 

Its components are to be added in the differential equation governing the motion of the 
particle. We are not interested in the solution of the homogeneous equation (Eq. 12) since 
it would be periodic and would average out. 

Before proceeding further, let us try to derive the transformation rule for vectors between 
the reference frames. In the warped ring having the frame of reference X'Y'Z', 



h + ecos0*, 

h + ecos (k x x' + k y y' — ut) 



From this we derive, 

— = -ek x sin <p , 
ox 

— = -ek v sm4> , 

dz' _ 
~ ' 

A unit vector normal to surface described by Eq. ED is written as, 



^(dz'/dx'f + (dz'/dy'f + (dz'/dz'Y 



-dz 



dz' t dz' 



dx' ^ ^ dy' ^ dz' 



(15) 

(16) 
(17) 
(18) 



or, 



k! = i 



ek x sin <p* 



+ k 
From Fig. ^ 



1 + e 2 k 2 sin 2 6* 



1 



1 + e 2 k 2 sin 2 0* 



+ 3 



ek y sin <p* 



1 + e 2 k 2 sin 2 . 



(19) 



O'P = (x - x*)i + (y - y*)j + (z - z*)k 

= ix 1 + jy 1 + kz 1 , 

and 

\0'P\ = Jx\ + y\ + z\. 



(20) 
(21) 
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The unit vector is given by, 

O'P * a?i - ?/i f zi 



= 9. (22) 



These relations yield, 
xi efc x sin i 



\Ai + 2/i + zf 


\Jl + e 2 k 2 sin 2 0* 




e/c y sin 0* 


\/xl + yl + z\ 


^Jl + e 2 k 2 sin 2 0* 




1 



(23) 



(24) 



(25) 



xf + yf + zf J 1 + e 2 A; 2 sin 2 0* 
From eqns. (|23|) and (J25j) we obtain, 

— = — ek x sin , 

or, x\ = Z\ek x sin 0*. (26) 
Similarly, from ()24)1 and (|25jl . 

— = — e/c„ sin 0*, 

or, 2/1 = Zie/Cy sin 0*. (27) 
From eqns. (JUE1I2ZI), 

— 0* = fca, (2ie/c x sin 0*) + k y {z\ek y sin 0*) 
= Z\6 [k 2 sin 0* + ky sin 0*) 

= ZieA; 2 sin0*. (28) 



Following techniques used by Shu (J1984J), the amplitude in the radial direction can be 
assumed to have the form, a(x) = a$exp [i J k x dx}. To get the principal Fourier component 
of the undamped disk, we choose a(x) = aoexp [ik x x]. Identifying the wave amplitude with 
e, one can clearly write ao as ao = e (Q 2 — u 2 ). 

Therefore, from eqns. ()13|) . the force component for satellite forcing is given as, 

z sat = e (yi 2 — uj 2 ^J cos(k x x + k y y — cut). (29) 

To get the components of satellite forcing from eqn. (J2*9~|) in the frame of reference placed 
at the mid-plane of warped ring system, one has to use the properties given in eqn. 
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The components of the equation of motions of the test particle are then given by, 



— = -2fl-p + 3n 2 x-u 2 Xl - (n 2 -uj 2 )e 2 k x cos(f)sm(j)*/W J 



= 2fi— - v 2 Vl - (tt 2 - uj 2 ) e 2 k y cos sin <j>*/W, 

d 2 z 

— = -Vl 2 z - v 2 z x + (ft 2 - uj 2 ) e cos (j)/W. 





(30) 



(31) 



(32) 



Here, v 2 = AirGp is the vertical oscillation frequency due to the self-gravity and p is the 
mass density of the ring matter. The first terms in the R.H.S. of Eq. (30) and (31) are the 
Coriolis acceleration term. The second term of eq. (30) comes from the difference between 
the centrifugal acceleration of the particle and the centrifugal acceleration of the observer 
in rotating frame of reference. The vertical forcing is associated with the coordinates in the 
primed (warped) frame. We consider the case when eK <C 1. i.e., for wavelengths which 
are large compared to the amplitude of the bending wave. In this limit, W ~ 1 and = 
0* + zien 2 sin0*. This means cos0 ~ cos0* — zien 2 sin 2 0*. 

In Eqs. (I30H32J) . we considered a kinematic approach, and not the fluid dynamical ap- 
proach. That is, we are not using Navier-Stokes' equation. This is because in the C ring 
which is made up of the neutrals, particles collide with one another very scarcely, once or 
twice in each orbital revolution (This many not be true in the partially ionized rings like E, 
F and G, where charged particles coupled to the magnetic field gyrate and the gyro-radius 
is much smaller than the ring thickness reducing the mean free-time and mean free-path 
drastically). Thus we prefer to treat the C ring as a collisionless system with particles sat- 
isfying Boltzmann distribution. However, once we obtain the local shear, we can use this 
to compute the local dissipation when the mean collision time is supplied. This two-step 
process has been adopted in our study. 

We first re-rite the z-component of the equation as, 



where, a 2 = Q 2 + u 2 . It is easy to transform the equation in terms of the phase 0* by using 




(33) 



d/dt = -ud/d(j>* and d 2 /dt 2 = uj 2 d 2 /d(f)* 2 . This yields, 

d 2 z x Zl \ 2 e 2 (tl 2 - lu 2 ) k 2 ] e 2 (n 2 -LU 2 ) Zl k 2 cos 20* 



(34) 



<90* 2 uj 2 I 2 2u 2 

Close to the resonance orbit, Q ~ u and the parameter, 



e 2 k 2 (n 2 - uj 2 ) 



(35) 



could be treated as a perturbation parameter (Chakrabarti, (jl989|)). Defining, 
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V 



a 2 7 2 



Eq. ({31} becomes 
d 2 Zl 



S*2 



2 , 7 2 cos! 
-V z x H — 



To obtain a complete solution, we expand Z\ in powers of 7 2 as, 



^i = 4 0) + E7 24 4 ij . 



(36) 



(37) 



(38) 



i=0 



Upon substitution in Eq. EH and equating coefficients of the powers of 7 2 we obtain, up to 
powers of 7 4 , the solution of z as, 



e cos 4>* + h sin 7/0* 
sin 7/0* 



hi 1 
T6~ 



sin (2 + 77) 0* sin (2 — rj) 0* 



1+7/ 

sin (4 — 7]) 0* 



V 



+ 



(I+77) (1 - 77) 8(1-77) (2-77). 
From Eqs. E21and|I3 we can now obtain the solutions for x and y, 



h^ 4 sin (4 + 77) 0* 
~64~ [8(1+7/) (2 + 77) 

(39) 



x 



h sin 



+ 



v 2 ek v h 



2oj 2 



1 + 77(16(1-7/)) 
fi 2 /c^ 2 - (I-77) 2 



cos(l — 7/)0* 



l + 77(16(l + r/)) 



fi 2 /^ 2 - (1 + 7/) 2 

7 2 /16 + 77(512(2 + 7/)) 



COs(l + 7/)0* 



(l + 7/)[fi 2 /^ 2 -(3 + r/) 2 ] 
77(512(1 + t/)(2 + 7/)) 



cos(3 + 7/)0* + 



7 2 /16 + 77(512(2 -77)) 
(1 -T])[n 2 /uj 2 - (3-t/) 2 ] 

77(512(1 -77X2-77)) 



[fi 2 /cu 2 - (5 - r/) 2 ] 



cos(3 — 7/)0* 
cos(5 — 7/)0* 



and 

V = 



[W/lu 2 - (5 + r/) 2 ] 
00* \ 



2/i cos 



u 2 ek x hQ 



1 +7 2 /16(l + 7/) 



cos(5 + 7/)0* 
'1 +7 2 /16(l -7/) 1 

fi 2 /^ 2 - (1 - 7/) 2 (1 - 7/) 

sin(l + r/)0* 
1 



(40) 



sin(l — 7/)0* 



n 2 /oo 2 - (1 + 7/) 2 (1 + 7/) 

7 2 /16 + 77512(2 -7/) 

fi 2 /cJ 2 - (3 - 7/) 2 (l-7/)(3-7/) 

7 2 /16 + 77512(2 + 7/) 1 



2 /cj 2 - (3 + T]) 2 (l + 7/)(3 + 7/) 

77512(1 -77X2- 77) 1 



sin(3 — 7/)0* 
sin(3 + 7/)0* 



n 2 /uj 2 - (5 - t/) 2 (5-77X3-77) 

77512(1 + 7/)(2 + 77) 1 



sin(5 — 7/)0* 
sin(5 + 77)0* 



n 2 /iu 2 - (5 + T]) 2 (5 + 77)(3 — 77) 

cos(l - 77)0* cos(l + 7/)0* 7 2 cos(l + 77)0* 



v 2 ek y h 
2u 2 



(l-nf 

7 2 cos(3 + 77)0* 
16(1 + 77) (3 + 77)2 



(1 + 77)^ 



16(1 + 7/) 3 



+ • • • other terms of higher order 



(41) 
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Nume rical simulations shown by Chakrabarti ([19891 ) and Chakrabarti & Bhattacharyya 
(20011) tend to indicate that the shear due to variation of radial velocity component along 
vertical direction is significant. However, for the sake of completeness, the variation of vertical 
velocity component along the radial direction is also considered. Analytically, the magnitude 
of this shear S could be computed from, 

, 2 / < £ > \ 2 f<Z> s2 



where, we used only the significant components in a corotating frame and <> indicates 
that averaging is being done over a complete cycle <fi* . To show how the shear looks like 
analytically, we show here the nature of the first term of Eq. (23), 
< x > P 

— = Q' (43) 
where, 



+ 



+ 



hVtuj 



/2vrfi\ 
cos I I 



1 



u 2 ek x h 



\ UJ 

'l + 7 2 /16(l-r/) 



i-v) 2 + 



2 , 1 +'f/ 1G ( 1 + V), 1+v) 2 



%-(l 



vY 



(—77 sin27rr/) 
2tt(1 + 77)(1 -77) 
l + 7 2 /16(l + r/) 



l + 7 2 /16(l-r 7 ) , 1 
(— sin 2irr)) 



,.,2 



vY 



r/) 2 v ~ ' '" J 2tt(1 + 77)(1 - 77) 
7 2 /16 + 7 4 /512(2 -r/) 



;i-^){S-(3-^) 2 } 

7 2 /16 + 7 4 /512(2 + r/) 



(3- v y 



+ 



l + ^){S-(3 + ^) 2 } 
7 2 /16 + 7 4 /512(2 -r/) 



-77 sin 27T77) 



2^(3 + ^(3-7/) 



[S-(3-r/)2 
7 2 /16 + 7 4 /512(2 + 7/) 



(3-i7) s 



(l+r?) g-(3 + 77) 2 
7 4 /512(l-r / )(2-7/) 



(3 + r/)^ 



sin 27T77) 



27r(3 + 7/)(3-7/) 



(5-t/) 5 



2 [ -?7sin27r?7 
1 ^ \27r(5 + 77)(5 - 77) 



sin 2iTr] 



2vr(5 + 77X5-77) 



(44) 



and, 
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—huj 

Q = — — sin 2nr] 

2,71 



T] 2 sin 27X7] ( 2 — T] 2 



64 



2tt(1 — ?7 2 )(4 — ?7 2 ) \1 -7] 2 
(i] 2 + 8) sin27rr7 
^(1 — ?7 2 )(4 — /7 2 )(16 — 77 2 ) 



rfirf — 10) sin 2tvt] r) sin 27rr] 



87r(l-r] 2 )(A-r] 2 )(lQ-7] 2 ) 2tt(1 - r? 2 ) 



(45) 



During the excursion of the particle in epicyclic orbits as well as in vertical motion inside 
the ring, the particle is expected to collide with a neighboring particle at least once in each 
orbit and in this process the vertical component of momentum is transferred. In this context 
we wish to remind that it is not enough just to have a differential motio n, the collision must 



also takes place in order that the momentum is transferred. Shu (J1984J) already pointed out 
that an anomalous viscosity may be required if one has to use only the differential motion 
due to the Keplerian distribution. The contribution ip = — |fi to the r — component of the 
shear stress S r< j, = rjip = where rj is the co-efficient of viscosity. This contribution at 

the resonance radius that we are considering is ip = 4.2 x 10 -4 comparable to the shear we 
computed above using Eq. (23). However, not being a strongly collisional system, rj need not 
be isotropic, and thus is likely to be larger for vertical excursion (due to collision among disk 
particles inside the ring of few meters in height) than radial excursion. For an isotropically 
collisional system, the excursion in all directions is likely to be statistically identical. In the 
X — Y plane, such excursion due to the local temperature causes the epicyclic motion. For 
simplicity, we assume that the excursions in the vertical and horizontal directions are similar 
and the epicyclic radius is of the order of the half-height of the ring. Thus, the mean free 
distance in both the vertical and the horizontal distance is either the epicyclic radius or the 
size of the particle, whichever is bigger. In case of E and F rings (which are not considered 
here), the situation may change, and the Keplerian contribution could be more compared 
with the contribution from the vertical excursion. 

Our analytical solution may be used for the calculation of the wave profile as well as the 
damping length of the bending wave once we model the collisional process properly. Simon 
and Jenkins (1994), in the context of A ring obtained the relation between the coefficient 
of restitution and the optical depth. For the time being we assume that the collision is 
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h=1.6m 1.8 2.0 2.2 




2 3 
Surface Density (gm/cm 2 ) 

Figure 2. Variation of the damping length (Xd) as a function of the surface mass density for various ring heights (h). A 
horizontal line at Ar; = 85km is drawn for reference. Note that there is no solution at high surface density. The shaded box 
indicates a prescribed parameter zone for the C ring as obtained from our model. 

nearly-elastic, and later we shall show a posteriori that this assumption may be justified. 
For a wave profile, we first c omput e the amplitude of the wave at the location of the vertical 
resonance from Rosen et al f)l99lh . 

, , 334 



e TV 
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1/2 



m . 



(46) 



where, the surface density a is to be computed in c.g.s. units. Here r v is the location of the 
ve rtical resonance. The radial wavelength of the launched wave is computed from Rosen et 

209cr 



A, (r) 



km. 



(47) 



(r - r v ) 

To calculate the damping length, we use the following procedure: we advance the wave by 
5r, and compute the reduced energy density due to collisional transport at r 2 = ri + 5r 
from, 



E w (r a ) = E w (n) - ^-^cir, 
c g \ r l) 



(48) 
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Figure 3. The wave profiles are shown. The outer one (dashed) shows for parameters selected outside our prescribed shaded 
parameter zone as given in Fig. 1^1 and the number of complete waves is few hundred while the inner one (solid) showing distinct 
wave-profile is obtained from a set of parameters taken from the shaded parameter zone. Here, the number of complete waves 
is 28 as observed. 



where, the energy density of the wave is calculated using, 
E w {r) = l -e 2 m\r) 



(49) 



and the rate of dissipation of the wave is obtained from, 

E w (n) = u k (n) <S> 2 . (50) 

The instantaneous kinematic viscosity z/^(r) is given by, 

u k (r) = n(r)r max (R\R%) . (51) 

Here, R is the size of the largest particle and R ep i is the amplitude of the epicyclic motion 
which is sim ilar to the thickness of the ring i.e., h. The group velocity c g (r) is given by Shu 
et al Jiasd), 

TlGa 



oj(r) — mQ(r) 



(52) 
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where, u>(r) = Qp — fi(r), Qp = —2-n jT titan is the pattern frequency, same as the angular 
velocity of Titan. For -1:0 resonance, m = 1 was chosen. 

We repeat the above procedure, till the energy of the wave vanishes. The traversed 
length gives the damping length of the bending wave. Those parameters which yielded 
results similar to the observed damping length and wavelengths may be correspond to the 
actual parameters of the disk. 



3 RESULTS FOR VARIOUS RING CONDITIONS 

In this section, we report some results of our calculation for rings of different characteristics 
of the Titan -1:0 resonance location. Computation for other rings will be studied elsewhere. 
The free parameters we considered are the ring height h and the surface mass density a. 
In Fig. we show variation of damping length (km) as a function of the surface mass 
density a for various thickness (h, marked on each curve). The horizontal line is drawn to 
mark the observed damping length of 85 km (Rosen et al. 1988). This line clearly shows 
that a particular damping length may be observed for two different values of surface mass 
density. It also shows that the ring height cannot be more than 3.6m. The shaded region 
around a ~ 1.5 gm cm -2 includes those parameters for which the number of oscillations lie 
between 25 and 35 and the damping length between 80 and 90 km. 

In order to show that the higher of the two values of a is really the solution, in Fig. El we 
compare the wave profiles for the two surface mass densities such as o = 0.31 gm cm -2 and 
1.57 gm cm" 2 for ring- height of 2.3m for which the damping length is 85 km. It is clearly 
observed that the case of o = 0.31 gm cm -2 (dashed curve) is invalid because the total 
number of oscillations is 148, far too high compared to the observed number. On the other 
hand, the profile for o = 1.57 gm cm -2 has 28 waves which is equal to the observed number 
of oscillations as given in Rosen et al (J1988J). 

So far, we obtained the results assuming the collision to be nearly-elastic. In reality the 
co-efficient of restitution e c need not be unity. Following simple kinematics, we can obtain 
the value of e c in terms of ring parameters from: 



U k <S> 2 71 2 

1 - Wtw ■ (53) 

After putting our obtained values of fj. ~ 1.2 x 10 -5 , h ~ 230 cm, r ~ o ~ 1.57 gm/cm 2 , 
VI ~ 2.835 x 10~ 4 , < S >~ 0.148 x 10~ 4 we obtain e c ~ 0.999297 which is very close to unity. 
It is to be noted that Simon and Jenkins (1994), using 6:7 Lindblad resonance parameters 
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obtained the e c = 0.9 for an optical depth of 1.57 (See, Fig. 4 of that paper). Though there 
is some discrepancy with our result, what is important is that the collision is nearly-elastic 
and thus the ring parameters derived by us may be realistic. 



4 CONCLUSION 



Rosen et al (1988J) inferred from their numerical modeling that the persistence of the wave for 
over 28 cycles indicated that the scale height of the C -ring could not be more than 2 to 3m. 
Even the radio occultation fitting indicated it would be rather ^ 2.6m for damping length 
of the order of ~ 85 km. Our analytical calculation based on the shear component in the 
vertical direction shows a good agreement with the observed results. For a given height and 
damping length we observed that there were multiple solutions for the surface density a. We 
were able to eliminate one of the solutions by explicit computation of the wave profile. The 
solution with the lower surface density has far too many number of oscillations that what was 
observed. We find that the most probable value lies somewhere around a ~ 1.57 gm cm -2 
and height h ~ 2.3m. We th erefor e believe that the new component of shear which was first 
introduced by Chakrabarti (1989) for Mimas 5 : 3 resonance, could be operating in Titan 
— 1:0 resonance as well. In future, we plan to explore other rings and other resonances. 

It is to be cautioned that our results depends on the nature of the amplitude and wave- 
length at the launching site as predicted using thin ring approximation ( equations l4~6Tl4*Tj) . 
It is possible that when the effects of finite thickness is incorporated the parameters may 
vary slightly. Another limitation of our work is that we did not explicitly include dissipation 
in the equation of motion themselves to obtain the wave profiles. We solved the problem 
in a two-step manner: we first compute the shear assuming no dissipation, and then use 
that shear to compute dissipation a posteriori. Similarly, no effect due to finite size of the 
particles, such as collisional break-up, particle spin etc. were included. We plan to follow 
some of the works presented in the Introduction to incorporate these effects to remove these 
short-comings in future. 
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